Trace initial interaction from final state observable 
in relativistic heavy ion collisions 
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In order to trace the initial interaction in ultra-relativistic heavy ion collision in all azimuthal 
directions, two azimuthal multiplicity-correlation patterns — neighboring and fixed-to-arbitrary 
angular-bin correlation patterns — are suggested. From the simulation of Au + Au collisions 
at y^SNN = 200 GeV by using the Monte Carlo models RQMD with hadron re-scattering and 
AMPT with and without string melting, we observe that the correlation patterns change gradually 
from out-of-plane preferential one to in-plane preferential one when the centrality of collision shifts 
from central to peripheral, meanwhile the anisotropic collective flow 112 keeps positive in all cases. 
This regularity is found to be model and collision energy independent. The physics behind the 
two opposite trends of correlation patterns, in particular, the presence of out-of-plane correlation 
patterns at RHIC energy, are discussed. 

PACS numbers: 25.75.Ld, 13.85.Hd, 25.75.Gz 



The data from current relativistic heavy ion experi- 
ments show that a new form of matter — quark-gluon 
plasma (QGP) has been produced at RHIC UM- The 
anisotropic collective flow V2 is supposed to provide us 
the information at the early stage of collision. The suc- 
cessful hydrodynamic description Q on observed mass 
dependence of V2 at low transverse-momentum range 
(pt < 2GeV ) shows that the observed dense matter be- 
haves like a perfect fluid rather than an ideal gas, and is, 
therefore, referred to as sQGP, although hydrodynamics 
can still not well fit all the observed data in the range Q 

A complete information about initial interaction and 
evolution is very important for correctly understanding 
the properties of the formed matter. The anisotropic 
collective flow Q V2 is the second Fourier coeflicient of 
the transverse-momentum distribution of final state par- 
ticles. It only globally characterizes the direction and 
strength of anisotropic distribution. The intrinsic inter- 
action (correlation) of final state particles is absent in the 
measure. 

Another related measure at the market is the 2-particlc 
azimuthal correlation @, 0] ■ It concerns the average cor- 
relation of two particles separated by a certain angle, no 
matter where the two particles are in the whole azimuthal 
space. The 3-particle, or 4-particle, azimuthal correla- 
tions focus on the same kind of correlations. How the 
particles in different azimuthal directions interact with 
each other can not be drawn from them either. 

Newly suggested two spatial-dependent correlation 
patterns, neighboring and fixed-to-arbitrary bin correla- 
tion patterns [1], provide information on how the parti- 
cles in different cells of phase space are correlated. They 
measure in the whole phase space the two typical corre- 
lations, i.e., correlations between local particles and be- 
tween particles with a certain distance in phase space. 
Therefore, they are good observable for tracing the ini- 
tial interaction among different azimuthal directions in 
relativistic heavy ion collision. 

In this letter, we will first introduce these two cor- 



relation patterns into relativistic heavy ion collision. 
Then, using the Monte Carlo models RQMD Q and 
AMPT [l^l, we demonstrate that the out-of-plane cor- 
relation pattern coexists with the in-plane one at RHIC 
energy, although the anisotropic elliptic flow V2 keeps 
positive there. Finally, the physical origin which gen- 
erates the two opposite correlations — out-of-plane and 
in-plane preferential correlation patterns is discussed. 

It is well-known that the general 2-bin correlation is 
defined as 
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1, 
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where mi and 7712 are the positions of the two bins in 
phase space and nm is the measured content in the mth 
bin. If there is no correlation between particles in the 
observed window, Cmi,m2 vanishes. 

We divide the 27r azimuthal angle equally into M bins 
and specify as the multiplicity in the mth angular bin. 



If we let TOi = m and 7712 = m + 1, C„ 



is reduced to 



the neighboring angular-bin correlation pattern, 

^m,m+l / 1 / 1 r. 
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If we fix the position of one bin, mi = mg, and vary the 
left one by m, it becomes the fixed-to- arbitrary angular- 
bin correlation pattern, 
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It is clear that these two correlation patterns measure at 
various azimuthal positions how the local particles and 
the particles separated by a certain angle correlate with 
each other. 

In order to apply these two measures to current rel- 
ativistic heavy collisions, we choose the RQMD and 
AMPT models as examples. As we know, the RQMD 
(relativistic quantum molecular dynamics) with re- 
scattering is a hadron-based transport model @. The 
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FIG. 1: The centrality dependence of neighboring angular-bin 
correlation patterns for Au + Au collisions at 200 GeV from 
the RQMD with re-scattering (first column), AMPT with 
(second column) and without (third column) string melting. 



final hadron interactions are implemented in the model 
by hadron re-scattering. Although the anisotropic col- 
lective flow produced by the model is much smaller than 
the observed data at RHIC, we can still see how the sug- 
gested correlation patterns behave for this kind of trans- 
port model. 

In contrary to the RQMD model, the AMPT is a 
"multi-phase" transport model, with hadron and parton 
interactions both included in. In the default AMPT, the 
transport in parton level is only partly, and is accompa- 
nied by string evolution. While in the AMPT with string 
melting, the parton level transport is fully taken into ac- 
count, and reproduces the observed anisotropic collective 
flow at RHIC So, using this model we will see how 
the azimuthal interaction changes when the anisotropic 
expansion develops from partly to fully. 

We generate 249,824, 204,004, and 686,278 events from 
RQMD and AMPT with and without string melting, 
respectively. The centrality dependence of neighboring 
angular-bin multiplicity correlation patterns for Au + Au 
collisions at ^snn — 200 GcV are shown in Fig. 1, where 
the results from the RQMD with re-scattering, from the 
AMPT with and without string melting are presented 
in the first to third columns, respectively. As an exam- 
ple, the fixed-to-arbitrary angular-bin multiplicity corre- 
lation patterns for Au-f Au collisions at the same energy 
from RQMD with re-scattering are shown in Fig. 2. Here 
we partition the whole azimuthal range 2tt uniformly 
into 50 equal size angular bins and in the first to third 
columns of Fig. 2 the positions of the fixed bins are lo- 
cated at mo = 1,6,12, corresponding to the azimuthal 
angles cf> = 0,j,^, respectively. The errors are statisti- 
cal only and most of them are smaller than the symbol 
size in these two and following figures. In this analysis, 
the (f) = refers to the direction of reaction plane in nu- 
clear collisions. In real experimental data analysis, it has 
to be determined event-by-event. 

We can see from Fig. 1 that the RQMD with re- 
scattering, the AMPT with and without string melting 
give qualitatively the same centrality dependence of az- 
imuthal correlation patterns. The correlation patterns 
are cos^ (j) like in peripheral collisions, then turn to flat 
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FIG. 2: The centrality dependence of fixed-to-arbitrary 
angular-bin correlation patterns for Au + Au collisions at 
200 GeV from the RQMD with re-scattering. 



in mid-central collisions, and become sin (j) like in near- 
central collisions. Here, we present only three centrality 
ranges to show their typical behavior. In fact, the corre- 
lation pattern changes gradually from cos^ (p like to sin^ (j) 
like with centrality. In the default AMPT the sin^ (f> like 
pattern is very weak at the centrality 40-60%, and be- 
comes more obvious in more peripheral collision with cen- 
trality 60 — 80% (not shown here). It is clear that two 
opposite trends dominate in peripheral and near-central 
collisions, respectively. In the mid-central collisions, the 
two trends turn to balance and the correlations become 
equal in all directions. Moreover, these important charac- 
teristics are independent of the specific assumptions im- 
plemented in the three models, in particular independent 
of the hadronization schemes assumed in these models. 

It is interesting and worthwhile to study these two op- 
posite trends in detail. In peripheral collisions, corre- 
sponding to the centrality at 40 ^ 60%, both correlation 
patterns show cos^ cj) - like patterns. The strengthes of 
the correlations are the largest at (j) = 0, tt, 2tt in neigh- 
boring angular-bin correlation pattern and are varying 
in fixcd-to-arbitrary angular correlation pattern — the 
nearer the fixed angular bin to = or tt, the larger 
the correlation strength. It tells us that the final state 
particles in the directions of in-plane {x-z plane, with 
z-axis corresponding to the beam direction and x axis 
along the impact parameter as shown in Fig. 3(a)) have 
the strongest neighboring correlations, and have also the 
strongest correlation with the particles in other direc- 
tions. But the strength of the correlation decreases when 
the fixed angle goes away from the direction of in-plane 
as shown in the first to third columns of Fig. 2. These re- 
sults show that the behavior of particles in in-plane direc- 
tions has the strongest influence to the behavior of par- 
ticles in all other azimuthal directions, in contrast with 
the fixed-to-arbitrary angular-bin correlation patterns in 
hadron-hadron collisions, where the back-to-back corre- 
lations are always the strongest one [ll|. These char- 
acters are similar to that of in-plane flow resulted from 
anisotropic expansion, which makes expansion privileged, 
or correlation stronger, in the directions of in-plane. 

On the contrary, approaching to the near-central col- 
lision, corresponding to the centrality at ^ 20%, the 
two correlation patterns show sin^ tf>-\ike patterns. The 
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FIG. 3: (a) An angular bin in the overlap zone, which is 
approximated as an ellipse in coordinate space; (b) The az- 
imuthal distribution of overlap area at four centralities. 



strongest correlations arc located in both ^ and i.e. 
out-of-plane (y-z plane) directions, instead of and tt (in- 
plane directions), and the particles in the out-of-plane di- 
rections have also the strongest correlation with the par- 
ticles in other directions. It indicates that the behavior of 
particles in out-of-plane directions has the strongest in- 
fluence to the behavior of particles in all other azimuthal 
directions. These characters are opposite to the in-plane 
flow but are similar to the out-of-plane one, which 

has been observed at Bevalac and SIS energies [3 ■ Such 
an out-of-plane flow was explained as that, before the 
anisotropic expansion, the participant nucleons, which 
are compressed in the overlap zone, cannot escape in the 
reaction plane (x-z) direction due to the presence of the 
spectator nucleons (squeeze-out effect jl4l. [l5j). produc- 
ing out-of-plane flow. However, when Lorentz contrac- 
tion effect becomes significant at ultra-relativistic colli- 
sions, the spectators will leave very soon from participant 
zone [l^ . The initial squeeze-out effect is supposed to be 
hardly realizable at RHIC energies Therefore, the 

observed out-of-plane correlation pattern is not due to 
the initial squeeze-out effect. 

It is clear that the anisotropic expansion and the late 
hadronization are impossible to produce stronge cor- 
relations in out-of-plane directions. Only the initial 
source anisotropy is preferential in these directions. The 
anisotropy of the initial overlap interaction area at non- 
central collisions results in a larger initial number of par- 
ticipant nucleons in the out-of-plane directions, which in 
turn generates stronger interaction in these directions. 

Let us estimate the azimuthal distribution of partic- 
ipants in the overlap zone and see whether it has the 
same centrality dependence as the out-of-plane correla- 
tion patterns. If we neglect the difference of density in 
radial direction and approximate the overlap zone as an 
ellipse, the number of participants should be approxi- 
mately proportional to the geometrical area of the sector 
in a given angular bin d^ as shown in Fig. 3(a), i.e., 

dA^part ^ ^ 1 2 ^ .4^ 

d./) d./) 2^ 2(62sin2</,-ha2cos2</.)' ^ ' 



where A^'part is the number of participant nucleons, iS^ 
the overlap area in the angular bin dc/), p the radius of 
the overlap zone in the considered direction, a and h are 
the semi- major and semi- minor axes of the overlap ellipse 
and can be deduced from the impact parameter 6i,„ and 
the radius R of colliding nucleus as a = y/ — ^hu/^j 
b = R- hj2. 

In Fig. 3(b), the azimuthal distribution of participants 
at four centralities are presented. We can see from the fig- 
ure that the azimuthal distribution of participants indeed 
behaves similar to sin^ 0, mimic the neighboring angular 
bin correlation pattern for near-central collisions shown 
in Fig. 1. The largest number of participant nucleons 
are located at and The amplitude of the distribu- 
tion is small in peripheral collisions, increases with the 
increasing of centrality, and reaches its maximum ampli- 
tude at the centrality 20 — 40%. It goes down slowly with 
the further increasing of centrality. The amplitude of the 
distribution will vanish in head on collisions as there is 
no difference in all directions of overlap area in such col- 
lisions. We indeed observed that the correlation patterns 
tend to flat at very central (centrality at — 5%) colli- 
sions in the AMPT with string melting. Such centrality 
dependence of the azimuthal distribution of participants 
can provide us a good understanding in the centrality 
dependency of the out-of-plane correlation patterns. 

The observed correlation patterns present the interac- 
tions between the particles in different azimuthal direc- 
tions. They are resulted from the competition of all kinds 
of interactions after collision, in particular, the initial in- 
teraction and subsequent anisotropic expansion. As long 
as the system is not fully thermalized, the initial out-of- 
plane preferential interaction will compete with the sub- 
sequent anisotropic expansion. In peripheral collisions, 
the overlap zone is small and so is the number of partici- 
pant nucleons. but the difference between the minor and 
major axes of overlap ellipse is large, and so is the dif- 
ference of pressure gradient. In this case the anisotropic 
expansion dominates the final observable, and the effects 
of initial interaction in correlation patterns are hidden. 
In near-central collisions, the overlap zone becomes large 
and the difference between minor and major axes of el- 
lipse is small, so that the initial interactions are strong 
enough to show themselves up in final observable. This 
is why the out-of-plane correlation patterns appear at 
near-central collisions. 

Therefore, the azimuthal multiplicity-correlation pat- 
terns provide us the information not only on the sub- 
sequent anisotropic expansion but also on the full initial 
interactions, which are sensitive to the anisotropic shape 
of initial overlap zone. 

We can further see from Fig. 1 that (1) the values 
of out-of-plane correlation patterns are almost the same 
in different models; (2) the values of in-plane ones vary 
greatly from one model to another; (3) the oscillation 
amplitudes of both the two opposite correlation patterns 
are the biggest in the AMPT with string melting. The 
reason of (1) is that the anisotropy of initial interaction 
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FIG. 4: The centrality dependence of neighboring angular-bin 
correlation pattern for An + An collision at (a) 60 GeV and 
(b) 5 GeV from the RQMD with re-scattering. 



is mainly determined by the geometry of initial interac- 
tion zone, which is the same in three models. Meanwhile, 
the subsequent anisotropic expansion is realized by the 
late transportation, which is very different in three mod- 
els, and so arc the results of (2). In addition, at the 
same collision energy, the strength of initial interaction 
also depends on the interaction mechanism [l^] ■ The 
anisotropic expansion is built up on the initial out-of- 
plane preferential interaction. A strongc initial interac- 
tion with full transporting results in large correlations in 
both in-plane and out-of-plane directions, i.e., (3). 

In order to see the energy dependence of the correlation 
pattern, the neighboring angular-bin correlation patterns 
for Au -|- Au collisions produced by RQMD with hadron 
re-scattering at 60 GeV and 5 GeV are shown in Fig. 4(a) 
and (b) respectively. It is clear that the higher the colli- 
sion energy, the bigger the values and amplitudes of out- 
of-plane and in-plane correlation patterns. We can also 
see from the figure that the presence of out-of-planc and 
in-plane correlation patterns is independent of collision 
energy, but dependent on the centrality only. 

However, if the initial eccentricity in coordinate space 
has not been transformed into, or inherited to, the 
transverse-momentum of final state particles, two oppo- 



site correlation patterns will be absent in the final ob- 
servable. This is the case for RQMD without hadron 
re-scattering, where the transverse-momentum distribu- 
tion is isotropic and correlation patterns are flat for all 
ccntralitics. 

As we know, anisotropic collective flow, V2, keeps pos- 
itive in RQMD with re-scattering and in the AMPT 
with and without string melting at 200GcV, 60 GeV 
and 5 GeV, which arc in consistent with the correspond- 
ing data. The results from the suggested patterns show 
that even if we observe in-plane preferential transverse- 
momentum distribution, it does not mean that the cor- 
relations, or interactions, have the same preferential di- 
rections. 

To the summary, we suggest to apply the two 
azimuthal multiplicity-correlation patterns in ultra- 
relativistic heavy ion collisions. They are neigh- 
boring and fixed-to-arbitrary angular-bin multiplicity- 
correlation patterns. From the simulation of Au -I- Au 
collisions at 200, 60 and 5 GeV by using the MC mod- 
els RQMD with hadron re-scattering and AMPT with 
and without string melting, we observe, model and en- 
ergy independently, that the correlation patterns change 
gradually from out-of-plane preferential one to in-plane 
preferential one when the centrality of collision decreases 
from central to peripheral, while V2 keeps positive in all 
cases. The in-plane preferential correlation patterns are 
resulted by anisotropic expansion, and the out-of-plane 
ones are found to be likely caused by the initial interac- 
tion due to the eccentricity of initial overlap zone. The 
experimental observation of the out-of-planc correlation 
patterns predicted in the present letter will offer us a bet- 
ter understanding for the initial pre-thermalization inter- 
action and its competition with the subsequent collective 
expansion 0, [inland is, therefore, called for. 
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